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Simulation of Relativistic Jets and Associated Self-consistent Radiation 
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Plasma instabilities excited in coUisionless shocks are responsible for particle acceleration. We 
have investigated the particle acceleration and shock structure associated with an unmagnetized 
relativistic electron-positron jet propagating into an unmagnetized electron-positron plasma. Cold 
jet electrons are thermalized and slowed while the ambient electrons are swept up to create a par- 
tially developed hydrodynamic-like shock structure. In the leading shock, electron density increases 
by a factor of about 3.5 in the simulation frame. Strong electromagnetic fields are generated in 
the trailing shock and provide an emission site. These magnetic fields contribute to the electrons 
transverse deflection behind the shock. Our initial results of a jet-ambient interaction with anti- 
parallel magnetic fields show pile-up of magnetic fields at the colliding shock, which may lead to 
reconnection and associated particle acceleration. We will investigate the radiation in transient 
stage as a possible generation mechanism of precursors of prompt emission. In our simulations we 
calculate the radiation from electrons in the shock region. The detailed properties of this radiation 
are important for understanding the complex time evolution and spectral structure in gamma-ray 
bursts, relativistic jets, and supernova remnants. 



I. INTRODUCTION 

Particle-in-cell (PIC) simulations can shed light on 
the physical mechanism of particle acceleration that 



occurs in the complicated dynamics within relativis- 
tic shocks. Recent PIC simulations of relativistic 
electron-ion and electron-positron jets injected into an 
ambient plasma show that acceleration occurs within 
the downstream jet 0, Hiii, [13, [H [M HI ■ 
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In general, these simulations have confirmed that 
relativistic jets excite the Weibel instability, which 
generates current filaments and associated mag netic 
fields [H, nil and accelerates, electrons |l3l [28| . 
Therefore, the investigation of radiation resulting 
from accelerated particles (mainly electrons and 
positrons) in turbulent magnetic fields is essential for 
understanding radiation mechanisms and their ob- 
servable spectral properties. 

Recently, synthetic spectrum has been obtained us- 
ing RPIC simulations d [Iljlljp, HI in order to 
examine "jitter radiation" 

'WW 

Further investi- 
gations are required to understand radiation mecha- 
nisms for gamma-ray bursts and variabilities in radi- 
ation from AGN jets. 



II. RECONNECTIONS 

Recently, reconnection has been proposed for addi- 
tional particle acceleration mechanism for AGN jets 
and gamma-ray burst jets [iO, El, [13, H, • Var- 
ious reconnection simulations have been performed; 
RPIC simulations Jj, |7|. [3C)|. [39| , resistive relativistic 
MHD (RRMHD) HMl^Tand two-fiuid [M [13 
simulations. In addition, the Kelvin-Helmholtz insta- 
bility (KHI) may also lead to particle acceleration [l| . 

In order to investigate the evolution of ejecta and 
associated emission we inject jets containing a per- 
pendicular magnetic field and associated convective 
electric field {E — — x varying the magnetic 
field strength, i.e^ magnetization parameter a = 
/ {Annj^meC^) [1, [3l- These simulations are dif- 
ferent from the previous simulations where jets were 
injected into a perpendicularly magnetized ambient 
13[. We have investigated the evolution of 



plasma 

colliding magnetized shells and calculate radiation as 
has been done theoretically and for RMHD simula- 
tions in order to include self-consistent microscopic 
effects m. 



Recent RPIC Simulations with Magnetized 
Jets Colliding with Unmagnetized and 
Anti-Parallel Magnetized Ambient Plasmas 



at a; = 500 A. The electron thermal velocity of jet is 
■^j^th — 0.014c, where c = 1 is the speed of light. Radi- 
ating boundary conditions were used on the planes at 
X — Xmin & aJmax- Periodic boundary conditions were 
used on all transverse boundaries. The ambient and 
jet electron-ion plasma has mass ratio mi/mc — 16. 
The electron/positron thermal velocity in the ambient 
plasma is v° = 0.05c and the ion thermal velocity is 
= 0.022c where c = 1 is the speed of hght. 
As in previous papers [11], the "flat" (thick) jet flUs 
the computational domain in the transverse directions 
(infinite width). Thus, we are simulating a small sec- 
tion of a relativistic shock infinite in the transverse 
direction. 

Figure 1 shows snapshots of the shocks generated by 
a jet propagating into an ambient plasma at simula- 
tion time t = 14500;";,^ with magnetization parameter 
aj'^^/cUp^ — 1. Here the jet carries a By magnetic field 
component with convective electric field component 
Ez. Panels in the left column show a case with no 
magnetic field in the ambient plasma, see the dotted 
blue line in Fig Ic. Panels in the right column show a 
case with an anti-parallel magnetic field (—By) in the 
ambient plasma, see the dotted blue line in Fig. Id. 

The anti-parallel magnetic field in the ambient leads 
to dramatic evolution in the collision region as shown 
on the right column. The electron density piles up 
at the jet front (Fig. lb), negative strong By (Fig. 
Id) and positive strong E^ (Fig. If) are found and 
indicate the occurrence of reconnection. In the case 
of no ambient magnetic field, jet electrons and ions 
propagate through the collision region (Figs. Ig and 
li), as opposed to the anti-parallel magnetized ambi- 
ent that hinders jet particle propagation through the 
ambient (Fig. lb). In the relatively short simulation 
time, electrons are accelerated promptly and strongly. 
As shown in Fig. 1, the magnetic fields play an essen- 
tial role in particle acceleration and, of course, in the 
generation of radiation. In this proposal we will sys- 
tematically investigate the effects of magnetic fields in 
relativistic flow collisions including reconnection. 



III. ELECTRON-POSITRON JET AND 
SYNTHETIC RADIATION 



We have performed simulations using a system with 
(Lx,Ly,L^) = (2005A,65A,65A) and a total of - a 
few million particles (8 particles/cell/species for the 
ambient plasma) in the active grid zones In the 
simulations the electron skin depth, Acc = c/wpc = 
10. OA, where Wpe = (47re^no/mo)^^^ is the electron 
plasma frequency and the electron Debye length Ac 
is half of the grid size. Here the computational do- 
main is three times longer than in our previous simu- 
lations [11]. The electron number density of the jet is 
as same as the ambient electron density and 7 = 15. 
In this study the jets collide with the ambient plasmas 



Figures 2a & b show the averaged (in the y~z plane) 
jet (red), ambient (blue), and total (black) electron 
density and electromagnetic field energy divided by 
the total jet kinetic energy from one simulation [25i] . 
The maximum density in the forward shocked region is 
about five times the initial ambient density. The jet- 
particle density remains nearly constant up to near 
the jet front. Current filaments and strong electro- 
magnetic fields accompany growth of the Weibel in- 
stability in the trailing shock region. 

The synthetic spectra shown in Figure 2c are ob- 
tained for emission from jets with Lorentz factors of 
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FIG. 1: Snapshots of the coUision of a perpendicularly magnetized jet with an unmagnetized (left column), or anti- 
parallel magnetized (right column) ambient plasma at simulation time t = 1450ct;j7c^ (Choi, Min, & Nishikawa 2011). The 
different panels show averaged values of the electron density (a) and (b), the magnetic field (c) and (d), the electric field 
(e) and (f), the phase space of electrons (g) and (h), and the phase space of ions (i) and (j). Reconnection occurs for the 
case involving anti-parallel magnetic fields and is indicated by the positive Ey component in Fig. If. 



7 = 10, 20, 50, 100, 300 and 1000 with cold (thin 
lines) and warm (thick lines) electrons (Nishikawa at 
al. 2011a,b,c,d,e). The radiation from the jet elec- 
trons shows a Bremsstrahlung-like spectrum for the 
eleven cases. 



However, it should be noted that at higher fre- 
quency the spectral slopes in Figure Ic are less steep 
than for a Bremsstrahlung spectrum. This is due to 
the Lorentz factor spread of accelerated jet electrons 
and a resulting higher average Lorentz factor. Addi- 
tional spectral extension to high frequency is due to 
electron scattering in the magnetic fields generated by 
the Weibel instability [H]. 



IV. CONCLUDING REMARKS 



The recent simulations of colliding jets into am- 
bient plasma with ant-parallel magnetic fields show 
that drastic evolution of the shock with piled-up mag- 
netic field with possible reconnection. At the colliding 
shock electrons are accelerated strongly, which may 
generate strong radiation. We will investigate this in- 
teresting evolution further including synthetic spec- 
tra. 
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FIG. 2: Averaged values of (a): jet (red), ambient (blue), and total (black) electron density, and (b): electric (red) and 
magnetic (blue) field energy divided by the jet kinetic energy at t = 3250 lo'^. (c): the spectra are for jets injected with 
Lorentz factors of 7 = 10, 20, 50, 100, 300 and 1000 with cold (thin lines) and warm (thick lines) electrons. The low 
frequency slope is approximately 1. 
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